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ABSTRACT 

It appears inevitable that reionization processes would have produced large-scale 
temperature fluctuations in the intergalactic medium. Using toy temperature mod- 
els and detailed heating histories from cosmological simulations of He II reionization, 
we study the consequences of inhomogeneous heating for the Lya forest. The impact 
of temperature fluctuations in physically well-motivated models can be surprisingly 
subtle. In fact, we show that temperature fluctuations at the level predicted by our 
reionization simulations do not give rise to detectable signatures in the types of statis- 
tics that have been employed previously. However, because of the aliasing of small-scale 
density power to larger scale modes in the line-of-sight Lya forest power spectrum, 
earlier analyses were not sensitive to 3D modes with > 30 comoving Mpc wavelengths 
- scales where temperature fluctuations are likely to be relatively largest. The ongo- 
ing Baryon Oscillation Spectroscopic Survey (BOSS) aims to measure the 3D power 
spectrum of the Lya forest, Pp, from a large sample of quasars in order to avoid 
this aliasing. We find that physically motivated temperature models can alter Pp at 
an order unity level at k < 0.1 comoving Mpc" 1 , a magnitude that should be easily 
detectable with BOSS. Fluctuations in the intensity of the ultraviolet background can 
also alter Pp significantly. These signatures will make it possible for BOSS to study 
the thermal impact of He II reionization at 2 < z < 3 and to constrain models for the 
sources of the ionizing background. Future spectroscopic surveys could extend this 
measurement to even higher redshifts, potentially detecting the thermal imprint of 
hydrogen reionization. 

Key words: cosmology: theory - cosmology: large-scale structure - quasars: absorp- 
tion lines - intergalactic medium 



The temperature of the intergalactic medium (IGM) is 
largely determined by how and when the cosmic hydro- 
gen and helium were reionized. Measurements of the mean 
temperature of the IGM at redshifts 2 < z < 6 appear to 
be consistent with the paradigm that hydrogen was reion- 
ized at z ~ 10 and helium was doubly ionized at z ~ 3 



1 INTRODUCTION 



all of these studies agree on the trends nor on their inter- 
pretation. 




* mmcquinn@berkeley.edu 



In addition to raising the mean temperature of the IGM, 
reionization processes would have heated the intergalactic 
gas inhomogeneously. Once imprinted, an intergalactic tem- 
perat ure fluctuation would fade away over roughly a Hubble 
time jHui fc Haimanll2003l 'l. Inhomogeneities in the temper- 
ature of the IGM alter the Lya forest absorption because 
the ionization state of hydrogen depends on the temper- 
ature, because the gas distribution on < 100 kpc scales 
is smoothed by thermal pressure, and because of thermal 
broadening of the absorption features. However, the vast ma- 
jority of studies of the Lya forest have assumed that there is 
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an approximately power-law relationship between tempera- 
ture and density. Our study investigates the detectability of 
realistic models for temperature inhomogeneities. 

The level of intergalactic temperature fluctuations 
has been con s traine d by previous a n alyses of the forest. 
iTheuns et all <|2002h and iLidz et al.l (|2010h placed upper 
limits on the allowed level of these fluctuations in the 
context of toy inhomogeneous heating models. In addi- 
tion, the standard model for the Lya forest has been suc- 
cessful a t explaining the statis tical properties of th is ab- 
sorption (Miralda-Escudc ct al. 1996; Hcrnquist ct al. 1996; 
iKatz et all 1 19961 ; iDave et all 1 19991 iMcDonald et al.ll2005l ). 
This picture posits that the fluctuations in the 2 < z < 5 for- 
est were driven primarily by density inhomogeneities, that 
the gas temperature can be approximated as a power law in 
density, and that the photoionization rate was nearly spa- 
tially invariant. The success of this model implies that de- 
viations from it cannot be large. 

Several theoretical studies of the Lya forest have exam- 
ined the i mpact of t e mper ature inhomogeneities on various 
statistics. lLai et all (120061 ) investigated the effect of tem- 
perature fluctuations on the Lya forest line-of-sight power 
spectrum. Surprisingly, they found that models with 0(1) 
fluctuations in the temperature (AT/T ~ 1) that correlated 
over tens of comoving Mpc altered this statistic at only the 
few percent-level for k < 5 comoving Mpc" 1 compared to 
models with power-law temperature-density relations. 

Other statistics hav e been found to be more sensitive 
to suc h inh omogeneities . iMeiksinl (|2000t ) . ITheuns fc Zaroubil 
( 2000), and lZaldarriagal (|2002h advocated the use of wavelets 
filters as a means to spatially identify temperature inhomo- 
geneities, and they showed that wavelets are sensitive to 10s 
of comoving Mpc fluctuations with AT/T m 1 when applied 
to 10 high-resolution Lya forest spectra. iLee fc Spergell 
(2010) found that threshold clustering functions (a popular 
statistic in the material sciences) could discrimina t e betw een 
different allowed thermal states. IFang fc White] \ 20041) ar- 
gued that a three-point statistic which correlates the large- 
scale flux with the small-scale power could place a strong 
constraint on deviations from a power-law temperature- 
density relation. Lastly, IWhite et alj (|2010T ) showed that 
large temperature fluctuations could have a significant effect 
on two- and three-point functions estimated from correlating 
multiple Lya forest sightlines. 

Here, we quantify the impact of 10s of comoving Mpc 
temperature fluctuations on a diverse set of Lya forest 
statistics. We use both toy models as well as the He H 
reionization simulations presented in lMcQuinn et all (I2009T ) 
to understand these effects. Section [2] discusses the differ- 
ent nonstandard contributions to fluctuations in the forest, 
concentrating on inhomogeneities in the temperature. Sec- 
tion [3] describes the methods of analysis used in our study. 
Section [4] (and Appendix A) quantifies the impact of mod- 
els for temperature inhomogeneities on different Lya forest 
statistics. Section [5] discusses how fluctuations in the ioniz- 
ing background could also alter Lya forest statistics. 

We adopt a flat ACDM cosmological model con- 
sistent with the m ost recent cosmological constraints 
l|Komatsu et al.ll201ll ). and we henceforth will use "Mpc" 
as shorthand for "comoving Mpc." All of our Lya for- 
est calculations normalize to the mean flux, (F), values of 



iFaucher-Giguere et all (|2008l ). which measured (F) = 0.69 
at z = 3 and (F) = 0.39 at z = 4. 



2 BACKGROUND 

In the Sobolev approximation, the optical depth for a photon 
to be absorbed as it redshifts across the Lya resonance of 
hydrogen is 

~ i-i a, r 4 r_ 12 ^— j dv/dx , (i) 

where A(, is the gas density in units of the cosmic mean, T4 
is the temperature in units of 10 4 K, H(z) is the Hubble 
parameter, and dv/dx is the line-of-sight velocity gradient 
(which is equal to H(z)/(1 + z) in the absence of peculiar 
velocities). Equation ((1} assumes that the hydrogen is in 
photoionization equilibrium with a photoionization rate of 
r_i2, expresse d in units of 10~ 12 s~ 1 . Me asurements find 
T_i2 w 1 fe.g jFaucher-Giguere et~aLll200St ). The T 4 " - 7 fac- 
tor arises because of the temperature dependence of the H I 
fraction in photoionization equilibrium, and many of the ef- 
fects studied here derive from this dependence. 

The amount of Lya absorption is determined by inho- 
mogeneities in the gas density, the peculiar velocity field, the 
photoionization rate, and the gas temperature. The first two 
sources are the standard contributions that were included 
in most previous studies and that are primarily responsi- 
ble for the statistical properties of the Lya forest. Almost 
all Lya forest studies also have assumed that the tempera- 
ture follows a power-law relation in density parametrized as 
T(Ab) = To A 7- (or they have effectively assumed a power- 
law relation by using simulations with a spatially uniform 
ionizing background). This study focusses on the impact of 
temperature fluctuations around T(Aj,), and it also consid- 
ers fluctuations in the photoionization rate. In addition to 
modulating the amplitude of TL ya , inhomogeneous heating 
alters 10s of km s~ absorption features via its impact on 
the gas pressure and on the thermal widths of absorption 
lines. 

Large-scale temperature fluctuations should be present 
in the IGM as a relic of reionization processes. The reioniza- 
tion of hydrogen and helium were spatially inhomogeneous, 
with some regions ionized earlier and some later depend- 
ing on their proximity to the ionizing sources. Models for 
these processes predict order unity fluctuations in the ion- 
ization fraction of the species being reionized on ~ 10 Mpc 
scales and accom panying order unity fluctuations in th e 
temperature (e.g.. iTrac et al.ll200Sl ; iMcQuinn et all 120091 ). 
Temperature fluctuations should have been imprinted on 
the IGM during reionization because: (1) different regions 
in the IGM would have been ionized by different energy 
photons, and (2) a gas element's instantaneous temperature 
depends on its ionization history. Once a typical element 
was reionized and heated, it would have subsequently cooled 
adiabatically owing to the expansion of the Universe (with 
Compton cooling off of the CMB and line cooling playing 
a minor role). This cooling proceeded until the gas element 
reached the temperature at which heating from ionizations 
of the res idual bound electrons balanced the cooling from ex- 
pansion (|Miralda-Escude fc Reeslll994l ; iHui fc Gnedm]|l997l ; 
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iTheuns et al.l [2002) . At this point, the temperature differ- 
ence with initially cooler elements would have been erased. 
At z < 6, it takes on the order of a Hubble time to erase 
temperature fluctuations^ 

The Lya forest indicates that hydrogen was likely reion- 
ized at z > 6. For any physical radiation spectrum that 
reionized the hydrogen, the intergalactic helium would have 
been at least singly ionized simultaneously with the hy- 
drogen. At z < 4, the relic temperature fluctuations from 
z > 6 reionization processes would have largely faded away 
l|Hui fc Gnedinl 1 19971 ; iTrac et all l200Sf ). However, a hard 
source of ionizing photons is required to doubly ionize he- 
lium (normal stars cannot do it). The current paradigm 
is that the intergalactic helium was doubly ionized by the 
radiation from quasars at z ~ 3 (IDavidsen et al] 1 19961 : 
Schave et al] l200d: lAgafonova et all 120071: iMcQuinnl 120091; 
Becker et all l201ll; IShull et all l20ld ; iFurlanetto fc Dixonl 
20ld ; iDixon fc Furlanettol 120091 ). Models predict that at 
z ~ 3 He 11 reionization would have imprinted large inhomo- 
geneities in the temperature, much larger than those r emain- 
ing from hydrogen reionization (|Furlanetto fc Oh] [2008b; 
IMcQuinn etaHl2009h . 

2.1 Simulations of the IGM Temperature 

Using radiative transfer simulations, McQ uinn et all (I2009T ) 
made predictions for the level of intergalactic temperature 
fluctuations expected during He 11 reionization. These simu- 
lations are employed extensively in our new study. In these 
simulations, the first regions in which the He 11 was reionized 
were closest to quasars and, therefore, were ionized by softer 
photons and heated by less than 10 4 K. These regions subse- 
quently cooled with the expansion of the Universe. Whereas, 
the last regions to be ionized were more significantly heated 
by ionizations from a hardened radiation background. Thus, 
the temporal extent of this process resulte d in large fluctua- 
tions in the intergalactic temperature. See I Abel fc Haehneltl 
(|l999h for discussion of relevant radiative transfer effects. 

Figure Q] shows slices through three snapshots of a 
430 M pc He 11 reionization simulation from IMcQuinn et al.l 
i|2009l ). The left panels are the He in fraction, the middle 
panels are the temperature, and the right panels are the 
change in the Lya forest transmission relative to the case 
with the same mean relation between T and Aj,, but without 
temperature fluctuations. By the end of He 11 reionization 
in this simulation, the temperature at the cosmic mean den- 
sity fluctuated between 10 and 30 kiloK, resulting in « 10% 
fluctuations in the Lya forest transmission. 

We briefly summarize the characteristics of the 
temperature fluctuations in these simulations, but see 

1 To quantify this, w e have solved for the evolution of To as in 
iHui fc Gnedinl d 199711 after hydro gen and helium reionization, as- 
suming a flat specific intensity at the H I, He I, and He II ionization 
edges (but the results depend weakly on this assumption). If half 
of the gas were at T = 15, 000 K and half at T = 25, 000 K at 
z = 3, by z = 2 the amplitude of fluctuations in To would have 
been 70% of what they had been, half at 17, 000 K and half at 
12,000 K. If instead half the IGM were at T = 15,000 K and 
half at To = 25, 000 K at 2 = 4, the amplitude of temperature 
fluctuations at z = 2 would have been reduced to 40% of its initial 
value. 




0.5 1 10 17.5 25 -0.1 0.1 
xHelll T (kiloK) dTransmission 

Figur e 1. Slices through simulation Bl from IMcQuinn et al.1 
(2009). The left panels show the He in fraction, the middle panels 
the temperature, and the right panels deviations in the transmis- 
sion relative to the case without temperature fluctuations but 
with the same mean T(A b ). Each panel is 430 Mpc X 430 Mpc X 
8 Mpc in size, with (from top to bottom) volume-averaged He II 
fractions of xjjeii = 0.1, 0.5, and 0.9 at z = 5, 4, and 3.5, respec- 
tively. 
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Figure 2. Distri bution of gas elements i n the T-Aj plane from the 
Dl simulation in McQuinn et al a smaller box version of 

the simulation featured in Figurc[T] The thick solid curves arc the 
mean T(Aj,), and the contours enclose 33,67,90,99, and 99.9% 
of the grid cells. These contours enclose the regions that have the 
highest density of gas elements in log T - log Aj, . 
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Figure 3. Top Panel: Magnitude of temperature fluctuations, 
quantified in terms of the dimensionless power spectrum, AS, = 
fc 3 P-p/27r 2 , where Pt is the power spectrum of <5y = T/f — 1. 
The thick curves show Ay for different He III fractions in the 
fiducial simulation of He II reionization. The thin blue dotted 
curve labeled "homogeneous" assumes 7 m 1.6 and T « 10 4 K. 
Bottom Panel: A 2 = fc 3 Px /27r 2 , where Px is the power spectrum 
of the overdensity in X . The blue dotted curve is for the same 
temperature model as in the top panel, and the curve labeled 
"Toy T Model" is A 2 for 7 - 1 = 0.3, where half the IGM is 
filled with uncorrelated 30 Mpc spheres with To = 20, 000 K and 
the other regions have To = 10, 000 K. The dashed green curve 
labeled x e is A 2 for electron fraction fluctuations, and the curve 
labeled "Density" is this for density fluctuations. 



iMcQuinn et all l|2009l ) for additional details. The solid 
curves in Figure [2] show the mean temperature as a function 
of At at different times in a higher resolution version of the 
simulation in Figure Q] All of these curves can be approx- 
imated by a power-law temperature density-relation with 
(7 — 1} w 0.35. However, there is significant scatter around 
this relation, with the scatter increasing as He n reioniza- 
tion proceeds. Moreover, the scatter is larger at lower Aj,, a 
property that affects the observability of these fluctuations 
(SectiongJ. The contours in Figure [2] enclose 33, 67, 90, 99, 
and 99.9% of the grid cells. 

The top panel in Figure |3] plots the dimensionless 3D 
power spectrum of temperature fluctuations, calculated from 
different times during the He 11 reionization simulation. We 
use the notation A x = k 3 Px(k)/2ir 2 , where Px is the power 
spectrum of the overdensity in X. The power law- like scaling 
of A| at k > 1 Mpc -1 owes to the temperature fluctuations 
being highly correlated with the small-scale density fluctu- 
ations because of adiabatic cooling and heating. The shape 
of Ay at k < 1 Mpc -1 owes to the large-scale inhomoge- 
neous heating during He n reionization (top panel, Figure 
[3]). Fluctuations at the ~ 10% level are present even for the 
largest modes in our 430 Mpc box during the bulk of He 11 
reionization. 

The bottom panel in Figure [3] compares Ay from the 
snapshot with xueii = 0.5 (dotted curve) with other relevant 



sources of fluctuations. A common toy model for IGM tem- 
perature fluctuations takes half of the IGM to be at 10 kiloK 
and the other half to be at 20 kiloK. The curve labeled "Toy 
T Model" is such a model, in which 30 Mpc, 20 kiloK bubbles 
are placed randomly until they fill half of the volume. The 
ringing behavior owes to the artificial top-hat bubble mor- 
phology. The fluctuations in this model are a few times larger 
than in the simulation. Thus, a similar level of temperature 
fluctuations as in the simulations should be more difficult 
to de tect than in th i s mod el. This is an importa nt point be- 
cause |TheunselaD (l2002f l and lLidz et"aH i|2010t ) found they 
could rule out this toy model using a wavelet analysis (Sec- 
tion [XT) . The curve labeled "Density" is k 3 P A (k) /2n 2 at 
z — 3, where the non -linear density power-spe ctrum Pa(&) 
is calculated using the lPeacock fc Doddsl i|l996h fitting func- 
tion. The temperature fluctuations in the simulation can be 
as large as the density fluctuations at the smallest k cap- 
tured, which suggests that such temperature fluctuations 
would have a significant effect on the large-scale transmis- 
sion fluctuations. 

He 11 reionization also produces additional free elec- 
trons, increasing the electron abundance (and TL ya ) by 8% in 
He ill regions. The power spectrum of the electron fraction 
from the xueii = 0.5 snapshot is represented by the dashed 
green curve in the bottom panel in Figure [3] The temper- 
ature fluctuations in the simulation are almost an order of 
magnitude larger than the fluctuations in x e . 

An inhomogeneous ultraviolet background (which 
modulates r_i2) is the final potential nonstandard 
source of transmission fl uctuations i n the forest and i s 
also the most s t udied dZuol ll992aHbt iMeiksin fc Whit 
2004 ICroftl 12004 iMcDonald et all 120051: jFurlanettd |200~ 



Mesinger fc Furlanettoll2009l ; IWhite et alj|201oj ). The curve 
labeled "Intensity" represents an empirically motivated 
model for intensity fluctuations at 2 = 3 (see Section [5}. 
Thus, the intensity power can also be comparable to that in 
the density at the smallest k. 



3 METHODOLOGY 

We use two of the He 1 1 reionization simulations presented in 
iMcQ uinn et all i|2009l ). The radiative transfer calculation in 
these simulations was run on a 256 3 grid in post-processing 
using the density field from either a 190 or 430 Mpc cos- 
mological TV-body simulation. Both simulations employ the 
same model for the sources. In addition to the He 11 reion- 
ization simulations, we use three supplementary simulations, 
all initialized with different rand om seeds. We use a 4000 3 
iV-body simulation described in IWhite et all l|201fj| ). This 
simulation has dimensions of 750 Mpc/h and the dynamics 
were softened at 100 kpc/h to approximate pressure smooth- 
ing for ~ 10 4 K gas. The other tw o supplementa ry simula- 
ti ons are 25 Mpc/h Gadget-2 runs (Springcl 2005) described 
in iLidz et all (|2010t ) with either 512 3 or 1024 3 dark matter 
and SPH particles. The supplementary simulations are used 
to overcome the resolution and sample variance limitations 
of the He 11 reionization simulations, as well as to include 
hydrodynamical effects absent in the reionization simula- 
tions. All of these simulations were initialized with cosmo- 
lo gical parameters that a re consistent with the measurement 



logical parameters that ar 
of iKomatsu et al.l (|201ll V 
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We add fluctuations in the temperature to the supple- 
mentary simulations using both toy models and a method 
that employs the temperature information from the He II 
reionization simulations. The latter method relies on the 
property that locally the T-At relation is well-approximated 
by a power law even though globally a power law is a poor 
approximation. A local power law holds because neighbor- 
ing cells have roughly the same thermal history. Specifically, 
this method obtains To and 7 for a cell in the He 11 reion- 
ization simulation by fitting this power-law model to T and 
Ab in the surrounding 5 cells, and it does this along one 
long diagonal skewer that cycles through the He 11 reioniza- 
tion simulation box. It then uses the resulting To and 7 to 
transfer the temperature field to density and velocity skew- 
ers extracted from the supplementary simulations (which we 
make periodic with length equal to the box size) . 

Figure|4]examines the validity of this method. The three 
sets of 5 3 points with different markers represent the values 
of T and Ab for the cells in three randomly selected regions 
of size 3.6 Mpc in the fiducial He 11 reionization simulation. 
The lines are the best-fit power law to each set of points. We 
find that AT /T ~ 0.1 — 0.2 is typical of the error around 
the best-fit T — At relation, which is much smaller than 
the global dispersion in the temperature. Thus, the local 
power-law approximation appears to be justified. However, 
this method assumes that the correlations between the large- 
scale temperature and density modes are unimportant. This 
is motivated by the fact that the correlation between tem- 
perature and large-scale density is weak in the reionization 
simulations owing to the rare nature of the sources. We will 
elaborate on the applicability of this approximation later. 

We will often use our method to add the He 11 reioniza- 
tion simulation temperature field at one redshift to a snap- 
shot in the supplementary simulations from another redshift. 
This makes possible visual comparison of how the tempera- 
ture fluctuations affect the statistic in question throughout 
the reionization process without the confusion of cosmologi- 
cal evolution in the density field. In addition, there is uncer- 
tainty in exact ly when He 11 reio nization by quasars would 
have occurred: iLidz et all (|2010h argued from the inferred 
temperature his tory that the end was a t 2 « 3.5. Whereas, 
iMcQuinnl (|2009h and lShull et al.l (|201Ch argued that it was 
not complete until at least z = 2.7, citing the presence of 
He 11 Lya Gunn-Peterson troughs down to this redshift. 



4 STATISTICS 

4.1 Small-scale Line-of-Sight Power Spectrum 

The line-of-sight power spectrum of transmission fluctua- 
tions is the principle statistic that has been used to derive 
cosmological constraints from the Lya forest. Here we inves- 
tigate how temperature fluctuations affect this statistic at 
k ~ 0.1 s km -1 , and Section l4~3l examines larger scales. Stud- 
ies of the forest marginalize over To, 7, and (in some cases) 
the reionization redshift to derive cosmological constraints, 
but the impact of the thermal history has the potential to 
be more complicated. 

Thermal broadening is the most important 
temperature-dependent effect on the small-scale power 
spectrum. Pressure smoothing of the gas can also smooth 




Figure 4. Fits to the local T—A^, relation using the z = 3.0 
snapshot from the 190 Mpc He II reionization simulation. The 
lines represent the best-fit T— A(, relation in 3 separate, randomly 
selected regions of 5 3 cells (corresponding to 3.6 Mpc on a side). 
The markers represent the value of T and A;, in each set of 5 3 
cells. 



out the small-scale transmission fluctuations. However, 
pressure effects are not included self-consistently in any of 
the simulations analyzed here. There are physical arguments 
for why pressure smoothing is less important to include 
than thermal broadening: 

(i) iGnedin fc Hul | 19981 ) showed analytically that in an 
expanding universe the scales at which pressure damps the 
growth of linear modes are a couple times smaller than 
the Jeans wavelength and smaller than the scale at which 
therm al br oadening erase s fluct uations. Both lGnedin fc Huil 
jl99Sl) and lPeeples et al.l (|2010h confirmed with simulations 
that thermal broadening dominates the exponential damp- 
ing of the small-scale power in the forest. 

(ii) The timescale for a typical Lya forest absorber to re- 
lax to equilibri um after a heating e vent is comparable to the 
Hubble time |Gnedin fc Hullll998l L During He 11 reioniza- 
tion, the additional pressure smoothing from the associated 
heating would not have had a significant effect on the statis- 
tics of the forest if this heating occurred within a redshift 
interval of Az « 2 (see Fig. 29 in lLidz et al]|2010h . 

Heating will also induce velocit y gradients that bro aden the 
majority of absorption systems l|Theuns et al.|[2000t ) . The ef- 
fect of this broadening mechanism on the Lya forest power- 
spectrum has not been quantified in as much detail. This 
process is unlikely to lead to an exponential cutoff in the 
small-scale power like thermal broaden because systems with 
zero width in redshift-space after including this process can 
still exist. 

Thus, this study ignores the dynamical response of the 
gas, a common approximation in such analyzes. However, 
this may bias the interpretation in this section and that in 
Section [All towards favoring higher temperatures. A treat- 
ment that includes this effect would be challenging. 

The shape of small-scale Lya forest power spectrum is 
sensitive to the average temperature as well as the temper- 
ature distribution. If /1 of the volume is at temperature 
Ti and the rest is at temperature T2, thermal broadening 
results in the small-scale power spectrum achieving the lim- 



© 2009 RAS, MNRAS 000, 1 



6 M. McQuinn et al. 



2 



10" 



10" 



10"' 



10" 



-2 



El, 



10"' 



10"* 



10- 




Xj=0.99, <T>=1.8 MloK 
Xj=0.94, <T>=1.9kiloK 
Xj=0.76, <T>=1.7 kiloK 
xpO.56, <T>=1.6kiloK 
Xj=0.40, <T>=1.3 kiloK 
x.=0.18, <T>=1.0 kiloK 



k (s km" ) 



10" 



Figure 5. Top panel: The line-of-sight power spectrum of 8p at 
z = 3 calculated from the 2 X 512 3 particle 25 Mpc//i simulation 
and focussing on the small scale affected by thermal broadening. 
The D5 curve is from the fiducial He II reionization simulation, 
the other curves are power-law T — A b models, and the "aver- 
age of two" curve is the average of the 10 kiloK and 20 kiloK 
models with 7 = 1.3. Bottom panel: Pp B at 2 = 3 calculated 
using the temperature field from different simulation snapshots, 
labeled by their He III fraction (xi in the plot). In both pan- 
els, the black points w ith error bars are the measurement of 
McDonald et al ] j2000f ) and the orange are the measurement of 
CrofT^r*aTTl200*2T) . The temperature field is taken from the snap- 
shots at z = 2.7, 3, 3.3, 3.6, 4, and 4.5 in the simulation. The thin 
Gaussian curve in the top panel is the wavelet filter used in Ap- 
pendix A, arbitrarily normalized. 



iting form at high k of 

P l F° s {k) -> h exp (-a2\ k 2 ) + (1 - /i)exp (-aT 2 k 2 ) , (2) 

where a — 2kb/m p . In the limit that T2 > T\ and 0T2 k 2 3> 
1, the power spectrum is solely determined by the regions 
with T\. 

An additional consideration is that the small-scale 
power at z ~ 3 is dominated by slightly overdense regions 
because these are the regions that contribute the narrow- 
est absorption lines. At z — 3, simple tests show that power 
on the exponential tail of Pj? s is primarily from regions with 
Aj, ~ 2 — 3, and this characteristic density decreases with in- 
creasing redshift. The dispersion in the temperature result- 
ing from an inhomogeneous reionization process should be 
smaller in overdense regions than in underdense ones. This 
trend is evident in the simulations (Fig. [2]), and it acts to 
reduce the impact of temperature fluctuations on the small- 
scale shape of Pp B . 

Figures [5] and [6] plot the small-scale line-of-sight power 
spectrum of Sf for several temperature models at z = 3 
and 2 = 4, respectively. All curves were calculated from the 
25/h Mpc Gadget-2 simulations (Section [3 1, and the points 
with error bars are th e mea surements of McDonald et al.l 
(|2000h and ICroft et all (|2002l ). In the top panels, the green 
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Figure 6. Same as Figure \E\ but for z = 4 and calculated from 
the 2 X 1024 3 particle simulation. See the key in Figure [5] for the 
mean temperature of the reionization simulation model curves. 



dashed (blue dotted) curves are for an IGM with 7 = 1.3 and 
T = 10 kiloK (20 kiloK). If 7 = 1.3 and half of the IGM had 
T = 10 kiloK and the other half T = 20 kiloK, then the 
small-scale power spectrum would be given by the magenta 
dotted curve, which has a slightly flatter shape than the sin- 
gle T0-7 models. Th e the shape of th i s curv e appears to be 
inconsistent with the lMcDonald et~ al. (2000) data. However, 
a toy model for the temperature field that is better moti- 
vated by the He n reionization simulations is half of the IGM 
at To = 15 kiloK with 7 = 1.5 and half at T = 25 kiloK 
with 7 = 1.2. The difference in the small-scale power spec- 
trum between these two states is negligible. (Compare the 
teal dot-dashed and grey quadruple-dotted curves in the top 
panels.) Therefore, the fluctuations in this toy model would 
not be detectable in P]? s . 

The bottom panels in Figures and [5] are the power 
spectrum of models that use the temperature field of differ- 
ent snapshots from the 190 Mpc He n reionization simula- 
tion. The IGM is gradually heated throughout He n reioniza- 
tion, which results in the small-scale power decreasing with 
time, or equivalently, with He in fraction. The amount of 
evolution is sensitive to the simulation's initial conditions, 
especially with decreasing ShcIii- This simulation was ini- 
tialized with To — 10 kiloK and 7 = 1.3 at the starting 
redshift of z — 6. If the simulation were instead initialized 
with To = 15 kiloK, there would have been less evolution. 

At the end of He n reionization, the simulation temper- 
ature reaches values that are sl ightly hotter tha n sugg ested 
by the z = 3 mea surements of iMcDonald et al.l (|2000l ) and 
ICroft et al l (|2002i ). The x H eIII = 0.94 curve in Figure cor- 
responds to z = 3 in t he simulation. The measurement of 
IMcDonald et al l f2000) at z — 4 in Figure [6] appears to be 
more consistent with slightly cooler temperatures than their 
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Figure 7. Top panel: The PDF of the normalized flux. The blue 
dotted curve represents a model with 7 = 1.3, and the teal dot- 
dashed curve for 7 = 1.0. The red solid curve represents a model 
with 7 = 1.3, but where half the IGM has T = 10 kiloK and the 
other half To = 20 kiloK. The green dashed curve is this PDF for 
the temperature field from the He II reionizatio n simulati o n. Th e 
points with errorbars are the measurement of I Kim et al. I d2007f) . 
Bottom Panel: The relative difference between the flux PDF in 
these temperature models and in the 7 = 1.3 model. 

2 = 3 measurement. For reference, the XHein = 0.4 output 
is at 2 — 4 in the reionization simulation. 

There are several systematics that may bias these inter- 
pretations. The observations are biased in the direction of 
having extra power by instrumental noise and especially by 
metal lines, and the simulations are also biased in the direc- 
tion of extra power as previously described because pressure 
smoothing is not properly captured and the simulation gas 
temp eratures (p r ior to our post-processing) are on the low 
side. iLidz et all l |201Ch estimated that pressure smoothing 
could alter the power at k ~ 0.1 s km" 1 by a maximum of 
20% (see their Fig. 29), but its impact should be significantly 
smaller for the He 11 reionization scenarios considered here 
where the h eatin g has been more recent than in this case in 
ILidz et all (^Old ). Uncertainties in the mean flux are also 
important for the comparison at 2 = 4 l|Lidz et al. I l2010l l. 

In conclusion, measurements of the small-scale Lya for- 
est power spectrum reveal that 15 < To < 25 kiloK at z = 3, 
with tentative e vidence for an de c rease in temperature from 
z = 3 to 4 in the lMcDonald et al.l (|2000l ) measurement. How- 
ever, finer determinations would be challenging because the 
systematics are at the level of the model differences. We also 
find that it would be difficult to detect physically motivated 
models for temperature fluctuations with this statistic. 

Appendix A considers the impact of the temperature 
models considered here on wavelet statistics. Wavelets are 
also a measure of the small-scale power in the forest, and 
are better suited than the small-scale behavior of P l p s for 
detecting temperature inhomogeneities. Despite this advan- 
tage, Appendix A reaches similar conclusions to this section. 

4.2 Flux Probability Distribution 

Another statistic that is often applied to the Lya forest is 
the probability distribution function of the normalized flux 
(FPDF). A distribution of temperatures at a given At, will 



broaden the FPDF relative to the case of a single T — At, 
relation. In fact, measurements of the FPDF have claimed 
to find too many hot voids at 2 < 2 < 3 compared to cos- 
mological hydrodynamic simulatio ns with standard, h omo- 
geneous thermal histories. This led iBolton et ail |2008) and 
IViel et all (|2009h to argue that an inverted T — At, relation 
(7 — 1 < 0) was needed to reconcile this discrepancy. He 11 
reionization by the observed population of quasars cannot 
generate an inverted relation l|McQuinn et al.l 120091 '). How- 
ever, perhaps an inverted relation is not required to explain 
the data, but may rather be accounted for b y additional dis- 
persio n in the temperature as conjectured in lMcQuinn et ail 
(2009). Dispersion would have made some of the voids hotter 
and, as a result, more transparent to Lya photons. 

The top panel in Figure \7\ plots the FPDF at z = 3. 
The point s with errorbars are the measured FPDF from 
iKim et all \200ti ). The curves are this statistic calculated 
from 1000 skewers through the 25/h Mpc simulation. The 
red solid curve represents a toy model in which 7 = 1.3 
and where half the IGM has To = 10 kiloK and the other 
half To = 20 kiloK. The green dashed curve was calculated 
using the temperature at the end of the He 11 reionization 
simulation (where XHeiii ~ 0.94), and corresponds to the 
simulation snapshot in which the dispersion in temperature 
is approximately maximized. The blue dotted curve repre- 
sents a model with 7 = 1.3 and with To equal to the average 
in the xuain = 0.94 snapshot. However, the dependence of 
this curve on To is weak. The bottom panel in Figure [7] plots 
the relative difference between the 7 = 1.3 model and the 
other temperature models featured in the top panel. 

We next estimate the magnitude by which tempera- 
ture fluctuations alter the FPDF. To proceed, we write the 
FPDF, V, as a Taylor expansion in St p T- 

V(F) « T(F) + A(5 Tp7 ) + B(5^ p7 ) + ..., 
A = §Fr-V(l-T), 

where V is the FPDF without temperature fluctuations, 
5t P 7 is the fluctuation in y -0 - 7 around the average T-Af, 
relation, F is flux for the average T-Aj, relation, (...) rep- 
resents an average at F (roughly corresponding to fixed 
Ab), and t = — log(.F). Both A and B are of the order 
of V(F), and we find A w 0.3 and B « -0.1 at F = 0.7. 
The leading order contribution from temperature fluctua- 
tions about the mean T — At relation is B {o% p7 ) . Plugging 
in {o~t P 7} ~ 0.2 2 , which is characteristic of the simulation 
temperature model at F ~ 0.7, we estimate that tempera- 
ture fluctuations should alter V(F) at the several percent 
level. 

The size of the residuals in the bottom panel of Figure 
[7] qualitatively agree with the above estimate that temper- 
ature fluctuations should produce a percent-level change in 
V(F). The predictions are not in quantitative agreement 
(and sometimes has the incorrect sign) because we have not 
included the normalization to a single mean flux in our ana- 
lytic expressions, which fixes the first moment of the FPDF. 
We conclude that a different effect is r equired to cre a te the 
« 20% suppression of the FPDF that IBolton et ail (|2008h 
finds at F w 0.7. Additional dispersion in the temperature 
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(even with a much larger amplitude than in our He II reion- 
ization simulations) is unable to explain the discrepancy be- 
tween the observed and simulated FPDFs. 

Changing 7 can have a larger effect on the FPDF than 
adding dispersion to the temperature because, in this case, 
the leading order term is A{5t P 7). (Compare the 7 = 1.3 
with the 7 = 1.0 curve in Figure [7]) Pe rforming a m o re de - 
tailed comparison than considered here. lBolton et all (|2008l ) 
found that an inverted r e lation with 7 w 0.5 provided the 
best fit to the lKim et~aT1 (|2007h FPDF measurement. How- 
ever, an inverted relation may not be the only explanation 
and would require an unknown heating mechanism. The 
FPDF is especia lly sensitive to the accuracy of continuum 
fitting (|Leell201lT ) and the efficacy of metal-line removal. 
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4.3 Large-Scale Correlations 

The remainder of this paper discusses large-scale correla- 
tions in the Lya forest. At scales where 5p = F/F — l and 8p 
are small, the 3D power spectrum of Sf can be expressed as 
bias factors times the power spectra of the different sources 
of fluctuations. In particular, 

iV(fc) « b 2 (c 2 PA{k) + 2Ge- 1 [P A T P 7(k) - P A j(k)} 

+e' 2 [PT P r(k)+Pj(k)]\ (4) 

where k is the Fourier wavevector, k — \k\, Tp7 is shorthand 
for the temperature field to the —0.7 power, J represents in- 
tensity, Pa x (fc) is the cross power spectrum between density 
and the overdensity in X, fi = k ■ h where n is the line-of- 
sight direction, and e « 2 - 0.7 (7 - 1). The G = (1 + gfi 2 ) 
factors arise from peculiar velocities (Kaiser 1987), and g re- 
flects that these distortions have a different bias than density 
fluctuations. Our simulations require g ~ 1 (Appendix B). 

The e and e 2 suppression of T- and J-dependent terms is 
approximate and results because the flux depends on At, T, 
and J via the combination Ap T~°' r J" 1 . This suppression 
holds in linear theory, but can be violated by higher order 
correlations. Equation Q suggests that at scales where ei- 
ther Pt P 7 or Pj are comparable to Pa, these fluctuations 
can have a large effect on Pp. 

4-3.1 Line-of-sight Correlations 

Previous Lya forest correlation analyses have primarily fo- 
cused on the line-of-sight power spectrum and neglected cor- 
relations between sightlines. The line-of-sight power spec- 
trum Pj? s can be expressed in terms of the 3D power spec- 
trum Pp as 

PF\h)= ^kPp(k,k [{ /k), (5) 

where fcy is a line-of-sight wavevector. Because fi = k\\/k, 
the projection to ID suppresses the impact of large-scale 
peculiar velocities (and the Pax terms; c.f. eqn.fi]). In ad- 
dition, 3D wavevectors with even k ^> ku still contribute to 
Pj? s at fcii . Figure[8]illustrates this aliasing effect, plotting at 
several fey the fractional contribution per log k to P]? s from 
different 3D modes with wavevector k. The red thick dashed 
curve represents this for fen = 10 _1 Mpc -1 and g — (the 



Figure 8. Fractional contribution per logfc to Pp B from different 
3D wavevectors with modulus. The curves (in order of decreasing 
dash length) are fcy = 10 -2 , 10 _1 , and 1 Mpc" 1 . The thick 
curves are for g = (no peculiar velocities) and the thin are for 
9 = 1- 
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Figure 9. The large-scale line-of-sight dimensionless power spec- 
trum for different temperature models at 2 = 2.8 (top panel), 
and the relative difference between these models and the 7 = 1.3 
model (bottom panel). The "Toy Model" curve assumes that 
7 = 1.3 locally and that 30 Mpc spheres with To = 20 kiloK 
fill half of the volume and the rest is at To = 10 kiloK. The er- 
ror ba rs are the SDSS measurement presented in McDonald et al. 
1120051 1. 

red thin dashed this for g = 1), approximately the value of 
the 3D wavevector where Pr is maximized relative to Pa in 
some of our models (Fig. EJ). Thus, Pp B (k n ) receives a sig- 
nificant contribution from approximately two decades in k. 
This aliasing dilutes the impact of large-scale temperature 
fluctuations! 2 ] 

Figure [9] quantifies the impact of temperature fluc- 
tuations on the line-of-sight power spectrum. The curves 

2 The calculations in Figures \E\ assume the simple form 
P F (k,k n ) = b 2 [1 + g(k }l /k) 2 ] 2 Pl n (k) exp[-fe 6 Tifc|/mj,], where 

Pj^ n is the linear theory overdensity power. We take g = 1 and 
T = 20, 000 K. We find that this functional form provides a decent 
approximation to the spectrum of Pp in our simulations. 
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are Pp, calculated using 22,500 skewers of side- length 
750 Mpc/h from the z — 2.8 snapshot of Run 1 in 
IWhite et al.l (2010) and using different temperature models. 
Each curve represents a different model for the temperature. 
However, the values of Pj? s in all of the considered models 
differ by less than 10%. These smal l differ ences qualitatively 
agree with the results of lLai et al.l (|2006T ). 

The points wi th errorbars in Figure [9] are the SDSS 
measurement from iMcDonald et al.l (|2005T ). The differences 
bet ween the plotted models are comparable or smaller than 
the IMcDonald et all (|2005l ) errorbars. Thus, the impact of 
these temperature models almost certainly could not be de- 
tected in P}? s once one marginalizes over the cosmology. 
Temperature fluctuations could bias cosmological parame- 
ter determinations from P]? s , but at no more than the few 
percent-level in the considered models. 

The calculations in Figure [5] implicitly assume that 
the large-scale temperature fluctuations do not correlate 
strongly with the large-scale density modes due to our 
method for adding temperature fluctuations. This is not al- 
ways a good approximation for the calculation that uses the 
He II reionization simulation temperature field. Accounting 
for these correlations would increase the difference with the 
D5 curve by approximately a factor of 2, as can be inferred 
from the discussion in the Section \4. 3. 21 

The next section discusses the 3D power spectrum of 
the flux, a statistic which avoids the aliasing issues in Pj? s . 
We show that temperature models that have almost no effect 
on Pj? s can alter the 3D power spectrum of the flux at an 
order unity level. 



4.3.2 3D Correlations 

The Sloan Digital Sky Survey Hi's Baryon Oscillation Spec- 
troscopic Survey (BOSS) aims to observe the Lya forest in 
the spectra of 1.6 x 10 5 quasars over an area of 8000 deg 2 . 
A future survey, BigBOSS, will push a magnitude fainter, 
increasing this sample by a factor of a few|j BOSS's high 
density of quasars will allow the first measurement of the 3D 
Lya forest power spectrum. A major goal of this effort is to 
detect the baryon acoustic oscillation (BAO) features. The 
3D power spectrum will also be useful for studying temper- 
ature and other nonstandard sources of fluctuations in the 
forest because it avoids the aliasing issues that diminish the 
impact of large-scale correlations in Pj? s . 

Figure \TU\ quantifies the effect of the temperature fluc- 
tuations in the He 11 reionization simulations on Pf- The 
thick set of curves are calculated from the 190 Mpc simula- 
tion, and the thin (shown in three of the panels) are from the 
430 Mpc one. The different normalization of the curves in 
the two simulations results primarily because of the different 
redshifts of the snapshots. The redshift of the 430 Mpc sim- 
ulation curves is Az ~ 0.3 higher than the quoted redshift 
because He 11 reionization occurs slightly earlier in this sim- 
ulation (owing to fewer recombinations in this simulation). 
When calculating these curves, we normalize to the observed 
mean flux at the respective redshift. The solid curves are Af. 



3 http://cosmology.lbl.gov/BQSS/ 
http : //bigboss . lbl . gov/lndex . html 
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Figure 10. 3D Lyo forest dimensionless power spectrum. The 
dotted curves are calculated assuming 7 = 1.3, and the solid 
curves adopt the temperature field predictions of the two He II 
reionization simulations. The thick curves are from the 190 Mpc 
simulation and are at the quoted redshift, and the thin curves 
(shown in three of the panels) are from the 430 Mpc simulation 
from a snapshot that is Az fa 0.3 higher. The vertical line seg- 
ments mark the centers of the three BAO features that upcoming 
quasar surveys are targeting as a standard ruler. The bottom pan- 
els also include a piecewise curve representing the measurement 
error multiplied by 5 of BOSS for the signal in each piecewise 
region as described in the text. The top panels include similar 
sensitivity curves but for BigBOSS. The cyan thin dashed curve 
in the bottom right panel also includes the prediction for Ai in 
a toy temperature model that has 30 Mpc spheres at twice the 
temperature filling half of the volume. 

computed using the temperature field from the He 11 reion- 
ization simulation at the quoted xuciu and, for comparison, 
the dotted curves show A F for a power-law T — At relation 
with 7 = 1.3. Changing 7 or To has little effect on the latter 
curve. 

Figure [TD] illustrates that temperature fluctuations can 
significantly alter A 2 F on large scales, changing its amplitude 
by as much as a factor of 2. In addition, the relative impact 
of temperature fluctuations at fixed rEHein largely agree be- 
tween the small- and large-box simulations. This agreement 
is suggestive that the impact of temperature fluctuations 
on these scales is not affected by how well the small-scale 
features in the forest are resolved. The bottom right panel 
shows for comparison A 2 F predicted for our toy model with 
30 Mpc bubbles at twice the temperature filling half of space 
in the larger box (cyan thin dashed curve). Interestingly, 
even though the temperature fluctuations are much larger 
in this model than in the simulations (see Fig.[3J), their effect 
on A 2 F is comparable to that in the simulation temperature 
models. 

The BOSS quasar survey is forecast to measure A 2 F at 
z rs 3 to ~ 10% accuracy in bins with Afc ~ k in the range 
10 -2 < k < 1 Mpc -1 , precision beyond what would be re- 
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quired to detect this extra power in some of our models. 
The bott om panels in Figure 1101 includes an estimate dis- 
cussed in iMcQuinn fe White! l|201ll ) for the measurement 
error of BOSS multiplied by 5, assuming 8000 deg -2 over a 
region of depth 500 Mpc. These curves are calculated as- 
suming n c ft = 3 x 10~ 4 Mpc" 2 for the z — 3 panel and 
h e g = 2 x 10~ 4 Mpc~ 2 for the z — 3.3 (corresponding to 
~ 2 — 3 per deg 2 ), where n e s is a noise we ighted number den- 
sity of quasars on the sky and denned in lMcQuinn fc White! 
(|201lh . The sensitivity scales as n~g . These number densi- 
ties take the quoted magnitude limit for BOSS of m = 22, 
assuming that all quasars below this magnitude are targeted 
and that S/N — 1 in 1 A pixels at m = 22. With these as- 
sumptions, BOSS will be several times more sensitive to Pf 
at z = 2 — 2.5 compared to at z = 3. In addition, the verti- 
cal line segments in this figure mark the centers of the three 
BAO features that the BOSS and other upcoming quasar 
surveys are targeting as a standard ruler. Temperature fluc- 
tuations have the potential to bias the measurement of this 
feature unless a fairly general functional shape is assumed 
for the continuum on which these features sit. 

BigBOSS aims to obtain a denser sample of quasars, 
and we forecast that it will be able to constrain Ap to z ~ 4 
if its survey strategy is similar to that of BOSS, and possibly 
to higher redshifts if this strategy is optimized to go deeper 
|McQuinn fc WhitdlioTlh . The dot dashed curves in the top 
panels in Figure[3]show what BigBoSS can achieve, assuming 
that it obtains spectra for quasars that are a magnitude 
fainter than BOSS such that n e g = 5 x 10~ 4 Mpc~ 2 at 
2 = 2.6 and 1 x 10~ 4 Mpc~ 2 at z = 4.2. 

Peculiar velocities are not included in the calculation 
of the curves in Figure \W\ and the gas is assumed to be in 
the Hubble flow. Peculiar velocities reduce the importance 
of temperature fluctuations. However, peculiar velocities do 
not contribute to the power in modes transverse to the line of 
sight, and, thus, we expect that the differences seen in Figure 
[I0]should be representative of the effect on transverse modes 
and the monopole component of Ap. Observations will be 
most sensitive to the Ap monopole (Section [6]). 

Figure QT] demonstrates that equation Q provides a 
reasonable description of the impact of temperature fluctu- 
ations in the reionization simulations. This figure compares 
the predictions for Ap from equation Q (red dashed curve) 
with that of the full numerical result (black solid curve). The 
red dashed curve is constructed by adding each contribution 
that appears in equation Q separately, and the value of b 2 
in this equation was derived by taking the ratio of Pf to Pa- 
The blue dotted curve is Ap for a power-law T — At relation 
with 7 = 1.3, the green long dashed curve and thick cyan 
dot-dashed are respectively (b/e) 2 A^ p7 and 26 2 /eAy p7A . 
(The thin cyan curves are — 26 2 /e Ay p7A .) The field St P 7 
is taken to be the overdensity in temperature around the 
mean T — At relation, a slightly different definition than 
used previously and that is motivated in Appendix C. The 
sum of the blue dotted, green long dashed, and ± the cyan 
dot-dashed curve yields the red dashed curve, which is re- 
markably similar to the black solid curve. (See Appendix C 
for a complementary analytic description of the impact of 
temperature fluctuations.) 

The history of Pt P 7A in the He II reionization simula- 
tion is complicated. This function is less than zero at the 
beginning of He II reionization and at the largest-scales be- 
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Figure 11. Comparison of the predictions for Ap using equa- 
tion [4] (red dashed curves) and the full numerical result from 
the 190 Mpc box (black solid curves). The blue dotted curves 
are Ap for a power-law T — A;, relation with 7 = 1.3, and the 
green dashed curves and cyan dot-dashed curves are respectively 
the contribution proportional to Pt p 7 and \Pt p 7a\- In particular, 
the thin cyan dot-dashed curves are —2 b e Pt p 7A an d the thick 
are 2 fee - 1 Pt p 7A- In each panel, the sum of the blue, green, and 
± the cyan dot-dashed curve yields the red dashed curve. 

cause quasars are heating their immediate surroundings (top 
left panel, Fig. [TTJ. Later, it becomes positive as the large- 
scale void regions are ionized. Constraining this nontrivial 
evolution from observations of Ap would inform He 11 reion- 
ization models. 

The contribution from Pt P 7A to Ap is of comparable 
import to the contribution from Pt P 7 in the four cases con- 
sidered in Figure [TT] despite Pt P 7A being smaller than Pr P 7- 
This occurs because Pt P 7A is enhanced relative to Pt P 7 by 
a factor of 2 e in equation Q . This enhancement explains 
why the temperature fluctuations in the simulations produce 
a comparable change in Ap to that in our toy temperature 
model, despite the toy model having a much larger Pt- If the 
He 11 ionizers were more abundant than assumed here (as 
would be the case if the faint end of the qua sar luminosity 
function were steeper; iFurlanetto fc Ohll2008al ). Pt P 7A could 
be significantly larger, increasing the impact of temperature 
fluctuations on Pp. 

In conclusion, temperature fluctuations could have a 
significant effect on the 3D power spectrum in the forest 
and would be detectable in future 3D Lya forest analyses 
such as with BOSS. Our simple analytic description for their 
impact is extremely successful. 



5 INTENSITY FLUCTUATIONS 

Several studies have investigated the im pact of inte n- 
sity (J) fluctuations on the Lya forest dZud Il992allb| ; 
iMeiksin fc White! 120041 ; ICroftl |2004 iMcDonald et al.l 12009 ; 
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lFurlanetto|[2009l ; iMesinger fc Furlanettdl2009l ). Here we at- 
tempt to understand how inte nsity fluctuations c ould affect 
3D correlati ons in the forest. iMeiksin fc White! (j2004T ) and 
ICrofj |2004 ) concluded that J fluctuations have a small ef- 
fect at 10—100 km s _1 separations in the Lya forest correla- 
tion function, but found that they could change this statis- 
tic at the 10 s of pe rcent-level on larger scales. Interestingly, 
IWhite et all (2010) found that intensity fluctuations could 
have an order-unity impact on 3D correlations. A detec- 
tion of intensity fluctuations would constrain the rarity of 
the sources of ionizing photons, and it would constrain the 
contribution of quasars versus that of galaxies to the meta- 
galactic ionizing background. 

Intensity fluctuations (equivalent to F_i2 fluctuations) 
alter the 3D Lya forest power spectrum on linear scales as 

P F (k)xs P F {k) + b 2 [e~ 2 P J (k)-2e~ 1 (l + gLi 2 )P AJ {k)] , 

(6) 

(c.f. eqn.|4| where Pf(&) is the flux power-spectrum without 
J fluctuations. If we assume a Euclidean space in which all 
H I ionizing photons experience a single attenuation length 
A and treat the quasars as continuously shining lightbulbs, 
then 



Pj = 



arctan(Afc) \ ' 
Afc J 



(L 2 ) 
n(L) 2 



+ b 2 P A (k) 



(7) 



(which is similar to the expression for Pj in 

IMesinger fc Furlanettd feoO^ and 



Paj 



/arctan ( Xk ) 



Xk 



b q P A {k). 



(8) 



Here, Pj, is the power spectrum of intensity fluctuations, 
Paj is the cross correlation between intensity and density, 
b q and n are respectively the luminosity-weighted bias and 
3D number density of the sources, and (Afc) -1 arctan(Afc) 
is (up to a constant factor) the Fourier transform of 
r -2 exp[-r/A]Q 

Finite quasar lifetimes alter Pj when the lifetime t q 
is comparable to or shorter than A/c, or 300 Myr for A = 
100 proper Mpc. The variance of J does not depend on t q 
when normalizing to a single luminosity function, but fi- 
nite lifetimes can reduce large-scale J correlations because 
J becomes uncorrelated in regions separated by distances 
greater than > ct q . These effects can be included in equa- 
tion © by substituting a more complicated window func- 
tion for (Afc) -1 arctan(Afc). However, the primary effect is 
for Pj to become white at wavelengths greater than ct q 
rather than at ~ A. Thus, detecting intensity fluctuations 
on 100 Mpc scales would place interesting constraints on 



A commonly used formula for the correlation function of in- 
tensity fluctuations, £j, is 



f j(r) = 



(L> 2 r A 



with 



JO 



du ■ 



exp 



1 - 



(9) 



(10) 



sinh v 

and N = 4-7rA 3 ?i/3. Note that this formula for £j (derived by 
fairly complicated means in IZuolll992bl l is just the convolution 
of r -2 exp[— r/A] with itself (times a factor that depends on the 
luminosity function). Put another way, it is the Fourier transform 
of the Poissonian term in equation 10. 




k (Mpc -1 ) 




1(T 3 , 
1(T 3 



1(T 2 10" 1 
k (MpcT 1 ) 



10 L 



Figure 12. Top panel: Dimensionless power spectrum of intensity 
fluctuations to the 1/2 power. The thick dashed curves are the 
Poisson component, and the thin dashed are the clustered compo- 
nent. Both sets of curves assume A = 500, 300, 150, and 70 Mpc 
(from bottom to top), corresponding roughly to the mean free 
path at z = 2, 3, 4, and 5. For reference, the thick black curve is 
[A; 3 P a /2tt 2 ] 1 / 2 in linear theory at z = 3. The Poisson-component 
curves assume n = 10 — 4 Mpc -3 and (L 2 ) / (L) 2 30, which set 
these curves' normalizations. The clustering curves assume b q = 3 
and z = 3. Bottom panel: Same as the top panel but where the 
power from the proximity region is smoothed out for the Pois- 
sonian curves as described in the text. Note that 10 — 2 Mpc -1 
is approximately the smallest wavevector that can be estimated 
from the forest. 



quasar lifetimes. In addition, light travel effects alter the 
isotropy of the J fluctuations because the emission will have 
travelled further in the line-of-sight direction. This would 
result in additional /i-dependent terms in equation ©. 

The top panel in Figure [12] plots the Poisson and clus- 
tering components of Aj = [fc 3 Pj/2-7r 2 ] (the first and sec- 
ond terms in eqn. [TJl for infinite lifetimes. The thick dashed 
curves are the Poisson component for A = 500, 300, 150, and 
70 Mpc (in order of decreasing amplitude). For reference, 
these A corre spond roughly to its measu r ed value at z = 
2, 3, 4 and 5 |Faucher-Giguere et al.l 120081 ; iProchaska et al.l 
2009). Our calculations assume that the quasar luminosity 
function has the form $(L) ~ {[L/L,] a + [L/L,]' 3 )" 1 with 
a cutoff 3 decades above and below L t , where a = 1.5 and 
13 — 3.5. These choices result in {L 2 )/{L) 2 ~ 30. Changing 
both the upper and lower cutoff by a factor of 10 and in op- 
posite directions results in a factor of 3 change in (L 2 ) / (L) 2 . 
Furthermore, the curves in the top panel in Figure [12] use 
n — 10 -4 Mpc -3 , representative of quasars at z — 2 — 3. 
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The lHopkins et al l (|2007l ) quasar luminosity function yields 
n = {1.2 x 10" 4 , 9 x 10" 5 , 5 x 10" 5 , 2 x 10" 5 } Mpc -3 at 
2 = {2.5, 3, 4, 5}. 

For three of the four cases featured in the top panel in 
Figure 1121 the power in the Poisson component is always 
larger than the power in the density field (the black solid 
curve). However, the normalization of the Poissonian com- 
ponent of Aj is highly uncertain. The Poisson curves in the 
top panel of Figure [T5] ar e roughly a facto r of 3 below the 
correspond i ng est imate in lFurlanettol ((2009) in part because 
iFurlanettd (20091) used /3 = 2.9 for which (L 2 )/{L) 2 w 100 
with our luminosity cutoffs. Most investigations of J fluc- 
tuations on the forest have us ed values for j3 that are more 
similar to the value used her e l|Meiksin fc W hite 2004; Croftl 
l2004lMcDonald et alj2005h . However. IFurlanettd (120091) ar- 

;ued t hat /3 = 2.9 provides a better fit to the Hopkins et al.l 

20071) luminosity function. 

The thin dashed curves in the top panel in Figure [T3] 
represent the contribution to Aj from source clustering at 
z = 3 for b q — 3 and for the same A as the corresponding 
thick curve. These curves fall below the Poisson component 
with the same A. However, the clustering component of in- 
tensity fluctuations will likely be the dominant source of J 
fluctuations if galaxies are the source of ionizing photons (as 
several studies have argued might be the case at z > 4; e.g. 
iFaucher-Giguere et al.ll2008l ). Because the clustered compo- 
nent also affects the flux field via Paj, this will enhance its 
effect beyond its contribution to Pj and also allow this con- 
tribution to be separated via its distinct angular dependence 
(Section [5]). In addition, the power in the clustering compo- 
nent of J is always larger than that in the density — the 
standard source of fluctuations in the forest — at k < A -1 
since b q > 1. These scales become observable in the Lya 
forest at z > 3. 



5.1 Effect of Intensity Fluctuations on Pf 

The previous section showed that Pj can be comparable 
to Pa in the Lya forest. If our simple analytic expression 
is correct (eqn. |BJ), this implies that J fluctuations signif- 
icantly increase Pf beyond Pf- However, previous numer- 
ical studies of the impact of intensity fluctuations on P l p s 
have considered models where the power in J was compa- 
rable to or larger than the power in A;,, and these studies 



k > 0.1 Mpc -1 (e.g.. ICrofd 12001 ; iMcDonald et al.ll2005l) . 
Here we try try to resolve this apparent discrepancy. 

Figure \T3\ features numerical calculations of the impact 
of intensity fluctuations on A % . The thin black dashed curve 
is the contribution to the total power from intensity fluctua- 
tions using equation (O and the fiducial model for J fluctu- 
ations with A = 300 Mpc. To generate the J field, we have 
convolved a random distribution of quasars which have the 
said luminosity function with the function r~ 2 exp[— r/A]. 
This ignores the clustering contribution, which Figure [12] 
suggests is smaller but non-negligible in this model. The 
solid red curve is Af> without intensity fluctuations, and 
the dashed green is the fully numerical calculation of the 
effect of intensity fluctuations in this model. The blue dot- 
ted curve is the prediction of equation (0 - the addition of 
the dashed black and solid red curves. The analytic model 
does poorly at capturing the impact of intensity fluctuations 



0.01 



< 



0.001 



y= 1.3 
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b /e A, 
X = 100 Mpc 



0.1 1 

k(Mpc~') 

Figure 13. Impact of intensity fluctuations on A F at z = 3, 
calculated using the 190 Mpc simulation. The solid red curve 
represents this for a uniform intensity field, and the blue dot- 
ted represents this for a fluctuating field with A = 300 Mpc, 
<L 2 )/{L) 2 = 30, and n = 10 -4 Mpc -3 . The green dashed curves 
are the predictions by adding A 2 to Ap as specified by equation 
J6}. The dashed black curve is ft 2 /e 2 A 2 in this model, and the 
teal dot-dashed curve is the same as the green dashed but for a 
model with A = 100 Mpc. 



at k > 0.1 Mpc -1 , but does better at smaller wavevectors. 
The teal dot-dashed curve is the same as the green dashed 
but for a model with A = 100 Mpc such that the intensity 
fluctuations are larger and thus have a larger impact on A F . 

The poor performance of the analytic model at k > 
0.1 Mpc - owes to the diverging character of J in quasar 
proximity regions. Within the proximity region of a quasar 
(within the distance r p (L) = [An(L)/L n A] -1//2 ), the effect 
on A% will be suppressed beyond the prediction in equa- 
tion <(6j by the exponential transformation from TL ya to 
transmission. Because TL ya ~ (1 + [fp/r] 2 ) -1 , the small- 
scale transmission power from proximity regions is damped 
exponentially with a kernel similar to that of a Gaussian 
with s.d. w r -1 /y/2. The bottom panel in Figure [12] is the 
same as the top panel, but where for the Poisson term in Pj 
given by fc 3 ^ 2 arctan(Afc)/Afc has been forced to transition 
to a constant function of k with the transition occurring 
at k — r p (L) -1 / v / 2- This operation is meant to approxi- 
mate the convolution of arctan(Afc)/A& with the Gaussian- 
like shape of the proximity region. The power in J is sig- 
nificantly affected by this operation at the largest k that 
are shown. The k where the proximity-region damping oc- 
curs are roughly k > r p ((L 2 ) 1//2 ) -1 , which correspond to 
k > 0.1 Mpc -1 in our fiducial model. Smaller k correspond 
to where the analytic model provides a good description in 
Figure [H 

Density fluctuations become larger than those in J at 
k > 0.3 Mpc -1 in three of the Poissonian models in the 
bottom panel of Figure [T^] (and the intensity power is also 
suppressed by an additional factor of e 2 in A F )- Wavevec- 
tors with k > 0.3 Mpc -1 correspond roughly to those where 
Pj? s has been measured. As with temperature fluctuations, 
the impact of intensity fluctuations are further diluted by 
aliasing in Pj? s such that 3D fluctuations are better suited 
for detecting them. Thus, it is unlikely that the intensity 
fluctuations in three of these models would have produced a 
detectable imprint on Pj? s . In fact, we have calculated the 
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impact of our fiducial J model on P]? s in the 25 jh Mpc box 
and, because of this small-scale damping, found no appre- 
ciable effect owing to these damping and aliasing effects. 

At z > 4, it is conceivable that intensity fluctua- 
tions contribute a substantial portion of the power even 
in P X p s because both A and n are smaller. Interestingly, 
iMcDonald et alJ l|2005h constrained the z = 4 mean flux 
from power spectrum measurements to be 10 — 20% higher 
than t he more direct measurement in lFaucher-Gieuere et al.l 
(2008). Such a disparity could potentially occur if an un- 
accounted source of fluctuations were cont ributing to P^? s . 
Some of the intensity models considered in lMcDonald et all 
(2005) [their Fig. 5] and more generally in models where 
{L 2 ) /[{L} 2 X 2 n] > 10 2 Mpc are able to produce a tens of 
percent change in the power at all scales. 

There are several possibilities for how the impact of 
intensity fluctuations on Pf can be distinguished from the 
fluctuations in the temperature. Temperature fluctuations 
should peak in amplitude at z — 2.5 — 3 if He II reionization 
were ending around these redshifts. In contrast, intensity 
fluctuations are likely to have increased monotonically with 
redshift. In addition, the amplitude of Pj is set by the factor 
(L 2 ) /[{L} 2 X 2 7T,3d], which can potentially be constrained by 
ancillary observations. 



6 SEPARATING THE COMPONENTS OF P F 

The signal we are interested in has the form Pf — Po + 
2g (i 2 Pi + g 2 fj, 4 P2 in linear theory, where the Pi are functions 
of only k = |k| and are given by Po = b 2 Pa + b 2 e~ 2 [Pj + 
Pry?], Pi = 2 b 2 P A + 2 b 2 e' 1 [Paj+Paiv?], and P 2 = b 2 Pa. 
Thus, temperature and intensity fluctuations affect Po and, 
to the extent they correlate with density, Pi. (Light-travel 
effects can generate additional fj 2 " terms that we do not 
consider.) Here we investigate the sensitivity of a survey 
to the terms in this angular decomposition. Detecting these 
angular components may prove easier in the Lya forest than 
it has for galaxy surveys because the velocity bias g is likely 
to be larger in the Lya forest than in galaxy surveys (for 
which g = l/b g , where b g is the galaxy bias). 
The Fisher matrix for the parameters Pi is 



dP F {k) dP F (k) 



var[P F (fc)] dPi 



dPj 



(11) 



where the sum is over Mk independent elements in a shell in 
fc-space centered at distance |fe| from the origin. In addition, 
var[Pf(fe)] is the variance on an estimate of Pf in a pixel 
centered on k. 

Let us first approximate var[Pp(fe)] as isotropic, which 
is the relevant limit for BOSS and, more generally, for all 
surveys with less t han 100 Lya forest spec tra per degree 
j|White et al.ll2010l ; iMcQuinn fc White! 120111 ). In this limit, 
the error on a measurement of Pi is 



*(k) 



-1,1/2 



r[P F (fc)] 



1/2 



(1.5, 4.O3- 1 , 10.7 5 - 2 ). 

(12) 

Also informative is the cosmic variance limit in which 
var[P F (fc)] = 2 6 4 P A (fc) 2 (1 + g^ 2 f, and for g = 1 results in 



<r(k) 



: P A (fc) (2.1, 20.4,28.2). 



(13) 



In both limits, a survey detects Po at the highest signifi- 
cance. Pi reveals how the temperature and photoionization 
fluctuations correlate with those of density. Because BOSS 
can constrain Pf in bins of Ak ~ k at z ~ 2.5 at the 1% 
level, it should also be able to place 10%-level constraints 
on Pi. 



7 CONCLUSIONS 

Significant temperature fluctuations with AT /T m 1 were 
likely imprinted in the high-redshift IGM by reionization 
processes, and they would have lasted for roughly a Hubble 
time thereafter. Temperature fluctuations would have been 
imprinted at 2 ~ 3 if He n reionization were ending near 
this redshift and, thereby, affected the majority of Lya for- 
est spectra. In addition, if hydrogen reionization ended at 
z rs 6, remnant temperature fluctuations could still be ob- 
servable in the forest from redshifts as low as z — 4 (e.g., 
ICen et al.ll2009h . This paper investigated the detectability 
of these fluctuations. 

Previous studies have shown that the z — 2 — 4 IGM 
cannot have been half filled with ~ 10 Mpc p atches of tem- 
perature 20 kiloK and the rest at 10 kiloK l|Theuns et al.l 
l2002l ; lLidz et al.ll201oh . However, we showed that half with 
To = 25 kiloK and half with To = 15 kiloK is consistent 
with recent Lya forest analyzes. Such a temperature dis- 
tribution is also close to what is found near the end of 
He 11 reionizati o n in the radiative transfer simulations of 
iMcQ uinn et al. I (|200gT ). If fact, we showed that the tem- 
perature history in the fiducial simulation of He 11 reion- 
ization in Mc Quinn et al.l (|2009l ) is grossly consistent with 
previous Lya forest measurements, although the observa- 
tions tentatively favor less evolution in the mean temper- 
ature. The temperature fluctuations in this simulation do 
not alter the small-scale Lya forest power spectrum (Sec- 
tion |4T]), the PDF of the normalized flux (Section |4~2)) . the 
large-scale line-of-sight power spectrum (Section I4.3|l . the 
wavelet PD F (Appendix A), and the three-point statistic 
proposed in lZaldarriaga et alj|200il (Appendix A) in a man- 
ner that is distinguishable from a power-law T— At, relation. 
Thus, the temperature fluctuations produced during az~3 
He II reionization by quasars would probably have evaded 
previous searches and also likely would evade future searches 
using these standard statistics. 

Interestingly, we find that realistic models for large- 
scale temperature fluctuations could have a significant effect 
on the 3D Lya forest power spectrum. This statistic can be 
measured by cross-correlating multiple quasar sightlines. In 
the line-of-sight power spectrum, the aliasing of small-scale 
3D density power to larger scale line-of-sight modes dramati- 
cally suppresses the prominence of temperature fluctuations 
for physically motivated models. However, we found that 
physically motivated temperature models could impart an 
order unity increase in the 3D power at k ~ 0.1 Mpc -1 . 
We showed that the impact of temperature fluctuations on 
the 3D power spectrum could be understood with a simple 
analytic model. 

Intensity fluctuations could also alter the large-scale 
correlations in the Lya forest. These too would have been 
hidden in the line-of-sight power spectrum by aliasing ef- 
fects. At sufficiently large scales, intensity fluctuations will 
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be the dominant source of fluctuations in the forest, and 
the fluctuations will become more prominent with increas- 
ing redshift. The impact of intensity fluctuations on Pf is 
likely to be distinguishable spectrally from the impact of 
temperature inhomogeneities. 

BOSS is forecasted to measure the 3D flux power spec- 
trum Pf with percent-level accuracy at z ~ 2.5 in a fc-space 
shell with Ak ~ k. Such a measurement would constrain 
the effects discussed here, and we anticipate that these ex- 
tra sources of fluctuations can be distinguished from other 
uncertainties in Pf (such as in the cosmology) when their 
impact is larger than ~ 10%. In addition, BOSS is capable of 
placing interesting constraints on this statistic to redshifts as 
high as z ~ 3 and BigBOSS to z ~ 4. Deeper surveys on an 
8 — 10 m class t elescope can push this me asurement to even 
higher redshifts jMcQuinn fc WhitefeOllT ). Furthermore, P F 
can be separated into different angular components with 
BOSS and future quasar surveys, which will facilitate the 
separation of the density, intensity, and temperature contri- 
butions. The cross correlation of the Lya forest with another 
tracer of large-scale structure (such as a high-z galaxy sur- 
vey) could also enable this se paration (|Guha Sarkar et al.l 
l201ll ; iMcQuinn fc Whitell201ll ) . 

The impact of temperature and intensity fluctuations 
will complicate attempts to constrain cosmic distances us- 
ing the BAO features in the 3D Lya forest. Because there is 
not one unique template for the spectrum of T and J fluc- 
tuations, marginalizing over their potential impact would 
likely require a fairly general parameterization for the con- 
tinuum on which the BAO sits. Such a marginalization pro- 
cedure would reduce the sensitivity to cosmological param- 
eters. Temperature and intensity fluctuations may also im- 
part (via their correlation with density) a scale dependence 
to t he BAO amplitud e. 

ICen et all (|2009h found that temperature fluctuations 
from models of H I reionization had a larger effect on the 
line-of-sight Lya forest power spectrum than we have found 
in our temperature models. They examined P l p s at z = 4 
and z — 5, using simulations of hydrogen reionization. They 
found that temperature fluctuations resulted in an increase 
in P]? s by 5 - 10% at k ~ 10~ 2 s km" 1 and 20 - 30% 
at k ~ 10 -3 s km -1 - the latter being roughly an order 
of magnitude larger than what our temperatu re models 
produ ce at z = 3. The larger effect found in ICen et all 
(2009) likely owes to two reasons. First, if the hydrogen 
were reionized by num erous dwarf galaxies as is assumed 
111 the lCen et all |2009h simulations, there would have been 
a stronger anti-correlation after this process completed 
between T and the large-scale density at relevant scales 
than in our He 11 reionization simulations. The effect of 
temperature fluctuations on Pj? s could be significantly 
enhanced by such an anti-correlation, as can be noted 
from equation Q. Second, the density power decreases 
with increasing redshift, which results in the temperature 
fluctuations (which tend to have AT /T ~ 1 because of the 
photoionization physics) becoming relat ively large r . The 
impact of temperature fluctuations in the lCen et al.l (|2009T ) 
models should be even more dramatic on the 3D Lya forest 
power spectrum. 
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APPENDIX A: OTHER LINE-OF-SIGHT 
STATISTICS 

Al Wavelets 

Several studies have suggested using wavelet functions 
to search for temperature fluctuations in the Lya forest 
|Meiksinll2000l ; lTheuns et al.ll2002l ; rZaldarriagall2002l ). This 
method uses the wavelet property that they are localized 
in configuration in addition to Fourier space. The idea is 




Figure Al. PDF of the wavelet coefficients Al with L = 
500 km s — 1 for several temperature model s. The points wit h er- 
rorbars are from the 2 = 3 measurement of lLidz et alj ll2010h and 
the solid red curves represent the 2 = 3 snapshot in the simula- 
tion. 




Figure A2. Same as Figure [ATI but for z = 4.2. The teal curve in 
the bottom panel represents the 2 = 4 snapshot in the simulation. 



to convolve the Lya forest transmission field with a wavelet 
that is sensitive to the amount of small-scale power to search 
for spatial variations in the power and, thus, temperature 
fluctuations. 

We compare the wavelet predictions of dif ferent tem- 
peratu re models with the recent measurement bv lLidz et al.l 
This study applied this technique to 40 high reso- 
lution, high S/N Lya forest spectra spa nning 2 < z < 4.5. 
The wavelet filter that iLidz et all (|2010h used is 



ty(x) — C exp (ikox) exp 



x 
"271" 



(Al) 



ILidz et all l|2010h chose C such that EjL, A x 9(x) 2 /N = 
1 for Ax = 4.4 km s" 1 . ILidz et al.l (| 20101 ) found s„ = 
69.7 km s _1 and ko s n = 6 to be a good compromise be- 
tween maximizing sensitivity to temperature while mini- 
mizing the impact of instrumental noise. Conveniently, the 
Fourier transform of 9(x), which we denote as ^(fc), is 
a Gaussian centered around k — ko with standard devi- 
ation s„ . The full width half maximum of ^f(k) spans 
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0.06 < fc < 0.1 s km -1 , and ^(fc) is plotted in the top panels 
in Figures [5] and [6] (thin curves, arbitrarily normalized). The 
wavelet filter that was used in oth er wavelet s t udies of the 
forest is qualitatively similar to the lLidz et al.l (l2010h filter, 
and we expect our conclusions are robu st to the exact filte r 
choi ce (such as the cu rvature statistic of lBecker et al.|[201ll ). 
iLidz et ail (|2010l ) primarily analyzed the PDF of 



A L {x) = dx \<f{x) o S F {xf 



(A2) 



where " o " denotes a convolution, and the integral averages 
the convolved signal over L — 500 km s" 1 in order to reduce 
the noise. Thus, their Al was the average of the wavelet 
power in a ~ 10 Mpc region. The mean of this PDF is sen- 
sitive to the average temperature (because it is a measure 
of the average power within the wavelet bandpass), whereas 
the width is a measure of the spatial variance in the temper- 
ature. However, cosmic variance in the forest is the primary 
determinant of the width of this PDF, and temperature fluc- 
tuations would manifest as an excess in the width over what 
is expected from simple models for the IGM thermal state. 

Figures [All and [A2I plot the predictions for the wavelet 
PDF at z = 3 and z — 4.2 using different temperature mod- 
els. These curves were calculated from the same simulation 
skewers as the Pj? s curves in Figures [5] and [6] The points 
with error bars are the measured values from ILidz et all 
(2010). The top pan el in both figures e xplore toy power-law 
T - A b models. (See lLidz et~ai1 (2010) for a more extensive 
comparison of such models.) 

The wavelet PDF of the 7 = 1.3, T = 10 kiloK model 
is quite discrepant with that of the 7 = 1.3, 20 kiloK model 
(Fig. IA1[) . The data favors the latter of these two models, 
and would also be inconsistent with a 50% mix of both tem- 
peratures (which would be the average of these two curves). 
Section [2] suggested that a more realistic toy fluctuating- 
temperature model is half the volume with To = 15 kiloK 
and 7 = 1.5 and the other half with To = 25 kiloK 
and 7 = 1.2. Unfortunately, the top panel illustrates that 
To = 15 kiloK and 7 = 1.5 produces a very similar wavelet 
PDF to the T = 25 kiloK and 7 = 1.2 case. Thus, as we 
found for Pj? s in Section EU a 50% mix of these two mod- 
els would be difficult to distinguish from a single To and 7 
model. 

The red solid curves in Figures [A 1 1 and [A"2l represent the 
simulation result near the end of He 11 reionization. These 
curves have a very similar width to the other curves in their 
respective panel despite the fact that these curves include 
dispersion in the temperature. The bottom panels show the 
wavelet PDF at different times during the simulation. Inter- 
estingly, none of these PDFs are noticeably broader than the 
power-law T — Aj, case. Over the course of He 11 reionization 
in the simulation, the mean of the PDF shifts to smaller 
values owing to the heating of the simulated IGM. The teal 
and red curves represent the z = 3 and z = 4 outputs in 
the simulation. The mean of the PDF in the simulation ap- 
pears to evolve slightly more between these redshifts than 
the data, consistent with what we found in Section [4.11 

In conclusion, measurements of the wavelet PDF con- 
strain the mean temperature and place limits on tempera- 
ture fluctuations at the level of AT /T m 1 for slightly over- 
dense gas. However, temperatur e fluctuations at the s maller 
level found in the simulations of lMcQ et al.ll2009l do not 
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Figure A3. Three-point statistic outlined in Section IA2I with 
k\ = 0.1 and &2 = 0.2 s km -1 . The top panel shows this statistic 
for different power-law T — A;, models and for the temperature 
in the simulation snapshot with shcIII = 0.94. The bottom panel 
shows this statistic using the temperature from several snapshots 
of the He II reionization simulation. 



significantly effect the width of the wavelet PDF and would 
be extremely difficult to detect with this statistic. 



A2 Small-Scale Power 
Correlation 



Large-Scale Flux 



IZaldarriaga et all |200ll ) and iFang fc White] (|2004h dis- 
cussed a particular three-point statistic that they argued 
provides an excellent test of the paradigm that gravitational 
instability shapes the transmission fluctuations in the Lya 
forest. They advocated the statistic 



C(k) = 



Phs F (k) 



\P h {k)Ps F {k)] 1/2 '' 



(A3) 



where P$ F (k) is the power-spectrum of Sp, Ph{k) is this of 
h, and PhS F (k) is the cross power between h and Sf. Here, 
h(x) = S H (x) 2 and 5b = S F (k)W kuk2 (k), where W kl>k2 (k) 
is a band-pass filter that transmits at fci < k < k2 and 
tildes distinguish the Fourier dual. Thus, h is the square 
of the band pass-filtered flux field, making it a measure 
of the bandpass power. This statistic is normalized such 
that a perfect correla tion between h with Sf yields (7=1. 
IFang fc White] (|2004h showed that this statistic could place 
strong const raints on the level o f tem peratu re fluctuations. 

As in [Zaldarriaga et al.l (|200lh and IFang fc White] 

(|2004f ), we take W kl ,k 2 (k) to be unity for ki < k < k^ and 
zero otherwise. The statistic C(k) is a measure of the cor- 
relation between Sf and the small-scale power in the high 
pass- filtered Lya transmission field. If gravity dominates the 
fluctuations in the forest, there will be more structure in a 
large-scale overdense (opaque) region, driving C(k) negative. 
Predictions for C(k) are shown in the top panel of Figure 
IA3I for different power-law T — At models. The curves are 
calculated from the 1000, 25/h Mpc random skewers drawn 
from the 2 = 3 snapshot of the hydrodynamic simulation 
and with ki — 0.1 and k-z = 0.2 s km -1 . These curves all 
have —0.6 < C < —0.5 at fc < 0.01 km s _1 , and the curves 
with the largest 7 have slightly smaller \C\. This trend re- 
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suits because, the larger the deviation from an isothermal 
relation, the more decorrelated a la rge-scale density mode i s 
from the small-scale density power ijZaldarriaga et al1l200lf ) . 

Temperature fluctuations from reionization could also 
decorrelate the small-scale density power from the large- 
scale flux. However, temperature fluctuations could also en- 
hance the negative correlation because a large-scale hot re- 
gion has increased transmission and also less small-scale 
power. We find that the former effect is most important 
in our reionization simulations: These temperature fluctua- 
tions decrease |C|, but not significantly. The red solid curve 
in the top panel of Figure [A"3l uses the temperature fluctua- 
tions from the XHein = 0.94 snapshot, and we find that |C| 
is slightly smaller for this case compared to the power-law 
T — Ab models. The bottom panel in Figure [A3l plots C using 
the temperature field from different snapshots of the He n 
reion ization simula t ion. 

iFang fc White] | 2004) showed that adding a lognormal 
dispersion at the Jeans scale around the mean T — Ab re- 
lation with standard deviation 0.2 results in |C| becoming 
significantly smaller. They used this result to constrain the 
standard deviation in this relation to be < 0.2 from a sin- 
gle quasar sightline. The temperature fluctuations in our 
He n reionization simulations have a standard deviation of 
<tt ~ 0.1 (Fig. [2}, but the simulations' temperature fluctua- 
tions are correlated over 10s of Mpc: In a large-scale hot or 
cold region, the dispersion in T — Ab is much smaller (Fig. [4]). 
Likely because of these large-scale cor relations, we find a sig - 
nificantly smal l er sup pression than in IFang fc White! J2004h . 
IFang fc White! l|2004f ) found the most dramatic effect for a 
higher pass filter than used in this section, in particular with 
fci = 0.2 and fc2 = 0.3 s km -1 . We find that if we use this 
filter choice, \C\ is reduced by a larger factor than for our 
fiducial filter (to a value as low as —0.3) but never to zero. 
However, even the choice k\ — 0.1 and k2 = 0.2 s km -1 is 
pushing the limit on wha t can be applied to even the highest 
quality Lya forest data l|Lidz et al.ll2010l ). 
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Figure Bl. Ratio of the flux power spectrum for modes with 
< fi < 0.25 to those with 0.75 < fj, < 1.0. The curves are mea- 
sured from the 25/h Mpc simulation at z = 2,3 and 4, assuming 
7 = 1.3. The increase in these curves at small-scales owes to ther- 
mal broadening. The approximate plateau in the ratio provides a 
measure for g. The horizontal lines are the expected large-scale 
anisotropy for g = 0.5 and g = 1.0. A larger box size would be 
required for a more precise determination. 

linear theory prediction in part because the absorption sat- 
urates in regions. In the limit that all the absorpt ion is from 
satura ted lines with bias b, g — 1 /b. Interestingly, McDonald 
(|2003h measured g = 1.5 at z = 2, larger than the linear the- 
ory prediction of g ~ 0.5. ISlosar et al. (2009) found a value 
closer to unity from large box, low resolution simulations. 
We have also measured g from our 25/h Mpc, 2 x 512 3 hy- 
drodynamic simulation and find values that are consistent 
with g ~ 0.5 — 1 (Fig. IBlJl . although a larger box size is 
needed for a more precise determination. 



APPENDIX B: PECULIAR VELOCITIES 

This Appendix quantifies the impact of peculiar velocities 
on large-scale flux correlations. Peculiar velocities have two 
effects: (1) The large-scale peculiar velocity field results in 
a redshift-space compression so that more systems appear 
in regions with converging flows. (2) The nonlinear peculiar 
velocity compresses or dilates Jeans-scale dense regions in 
absorption space (and dilates Jeans-scale voids). On large 
scales, the latter effect enters by altering the bias of the for- 
est. The former effect produces the redshift-space anisotropy 
of this signal. It affects TL ya via the factor (H a + dv/dx) -1 , 
where v and x are the peculiar velocity and comoving dis- 
tance along the line of sight. In Fourier space this factor is 
simpler a nd approxima tely equal to 1 + fl m (z) 06 n 2 5, where 
(i = k ■ h (|Kaiserlll987l ). Thus, on linear scales the redshift- 
space power spectrum of 8f is likely to have the form 

P F « b 2 [G 2 Pa + 2 G e" 1 Pax + e~ 2 Px] , (Bl) 

where G = (1 + g^ 2 ), g is the large-scale bias of velocity 
fluctuations, and X is a placeholder for T or — J fluctua- 
tions. The average of G 2 over solid angle is 28/15. In linear 
theory, g « Q m (z) ' 6 e . However, g will depart from the 



APPENDIX C: A SECOND DERIVATION OF 
EFFECT OF INTENSITY FLUCTUATIONS ON 

Pf 

This section develops a more sophisticated understanding 
for how A% is altered by large-scale temperature fluctua- 
tions. Let us assume £t P 7(t, Ab) does not depend on Ab, 
where ^r P 7{r, Ab) is defined as the correlation function of 
T^ 0,7 /(T^ ' 7 ) — 1 and Ta is the mean temperature at a 
given Ab- This differs somewhat from our previous definition 
of the 5t P 7 field as the fluctuation field, but the distinction 
makes little difference in practice. 

We define to as the Gunn-Peterson optical depth at 
Ab = 1, and the value of to is set by the mean flux normal- 
ization. We also define X = A 2- ' 77 and the normalized flux 
as F = exp[— tqX(1+St)] - Finally, we neglect peculiar veloc- 
ities, and we assume that the temperature fluctuations are 
uncorrelated with those in density and that they are Gaus- 
sian. While we found in Section 14.3.21 that the correlations 
with density are important and it is likely that the tempera- 
ture fluctuations are not Gaussian, the results of this model 
are illustrative and generalize beyond these assumptions. 

With these assumptions and definitions, the correlation 
between the normalized flux in a region with density Ai and 
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a region separated by a distance r with A2 is 
(FxF^t = f d6 Tl dS T2 e -o(x l( i + * Tl ) + x 2 (i + ^ 2 )) 



l(5 Ti 5 T2 )C- 1 (5 Tl S T2 ) T 
X — (CI) 

v /(27r) 2 detC 



e -r„(X 1+ X 2 ) x exp[ Z0^pi (X 2 + X |)] 



xexp[r 2 A-iJf 2 ^ Tp 7(r))] , (C2) 
» FiA (1 + To 2 ^!^ ?T P 7(r)) , (C3) 

where F = exp[-T X], C = (1, ^Tpr/orpri £tW0t P 7i 1)j 
and represents an ensemble average with respect to 

y. To go from equation (|C2|) to equation (|C3|) . we used that 
£rp7(f) < 1 to expand the exponential. Furthermore, we 
set the second exponential term in equation (IC2[I to unity. 
Since regions that dominate the flux correlation function 
have tqX ~ f , the error from dropping this term is of order 
a T P 7 ^ 1- While this is comparable to the temperature term 
that we kept, it is less interesting because it does not depend 
on ^T P 7(r). In addition, most of its effect is re-absorbed in 
the renormalization to a single mean flux. 

To calculate the flux correlation function, we average 
(F\F2)t over Ai and A2, which yields 

i F = {F 1 F 2 ) TA ^^ F{r) +K 2 ^ T (r) + r 2 ^ F (r)^ T (r), (C4) 

where f p(r) is the unperturbed flux correlation function and 

K = t (F 1 X 1 )a 1 , (C5) 

= r [ dAA 2 - 7 ' 1 exp[-T A 2 - 77 ]p(A). (C6) 



The function p(A;,) is the volume-weighted density proba- 
bility distribution. We have dropped terms in equation (|C4[) 
that involve additional £(r) factors as well as the connected 
moments that are not incorporated in K 2 £,r(r). 

We are interested in the impact of temperature fluc- 
tuations on > 10 Mpc correlations. At these scales, it is 
justified to drop the last term in equation (|C4|) . With these 
simplifications, the power spectrum of fluctuations in the 
normalized flux is given by 

P F (k)^P p (k) + ^P T (k). (C7) 

Thus, the correction proportional to £t(t") has bias K/F 
in this model. We calculate K 2 / F 2 — 0.18 at z = 4 and 
0.04 at z = 3, assuming To = 2 x 10 4 K, 7 = 1-3, and 
r m = I0~ 12 s"\ using the lMiralda-Escude et ahl ||2000h fit- 
ting function for p(A). It turns out that K 2 /F 2 is almost 
identical numerically to the corresponding factor that ap- 
pears in equation Q, b 2 /e 2 . 



